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Abstract 



We study the equilibria of a self-gravitating scalar field in the region out- 
side a reflecting barrier. By introducing a potential difference between the 
barrier and infinity, we create a kink which cannot decay to a zero energy 
state. In the realm of small amplitude, the kink decays to a known static 
solution of the Einstein-Klein- Gordon equation. However, for larger kinks the 
static equilibria are degenerate, forming a system with two energy levels. The 
upper level is unstable and, under small perturbations, decays to the lower en- 
ergy stable equilibrium. Under large perturbations, the unstable upper level 
undergoes collapse to a black hole. The equilibrium of the system provides a 
remarkably simple and beautiful illustration of a turning point instability. 
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I. INTRODUCTION 



We describe here a simple model of a system exhibiting a gravitationally induced tran- 
sition from an unstable static equilibrium to a lower energy stable equilibrium (for small 
perturbations) or to a black hole (for large perturbations). The system is the general rel- 
ativists version of a spherically symmetric scalar wave $ reflecting off an inner mirror of 
radius R. By holding the mirror at a fixed potential with respect to infinity, we create kink 
boundary conditions which do not allow decay of the field to a zero energy state. In the 
nongravitating case, any initial configuration would decay to the static solution of the wave 
equation, $ = AR/r, where A is the potential of the mirror; and this static solution is 
stable under perturbations. In the self-gravitating case we will show there is a critical value 
of kink potential A c about which the system displays a classic example of a turning point 
instability [1,2], responsible for the above phenomena. Above this critical potential, there 
are no equilibria and all configurations collapse to a black hole. The behavior of this system 
illustrates the remarkable richness of the physical properties that result from the mixture of 
long range gravitational forces and nonlinearity. 

Our theoretical model is the null-timelike initial value problem for a spherically symmet- 
ric space-time satisfying the coupled Einstein-Klein-Gordon equations for a massless scalar 
field. We pose this mixed boundary value problem in the region of space-time outside a 
timelike inner boundary and to the future of an initial null hypersurface. The mathematical 
details of this problem have been described before in an investigation of the decay of scalar 
wave tails [3]. In that case, the scalar field was set to zero on the inner boundary, so that 
the boundary acted as a perfectly reflecting mirror. Here, by holding the inner boundary 
at a fixed potential, its role as a mirror is unchanged due to the local gauge invariance of 
the system with respect to the transformation <3> — > <3> + constant. However, the potential 
difference with respect to infinity introduces nontrivial global effects. 

We set the geometry of the mirror to correspond to a surface of radius R enclosing 
a flat region of spacetime. The evolution of the system then depends uniquely upon the 
characteristic data for the scalar field on the initial null hypersurface, which extends from 
the mirror to future null infinity T + . 

We describe this system in terms of a null coordinate system with the spherically sym- 
metric line element 




(1) 



r 



In these coordinates, the Einstein-Klein- Gordon equations reduce to [4,5] 




(2) 




(3) 



and the scalar wave equation □$ = 0, which takes the form 




(4) 
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The initial null data necessary for evolution consists of <&(tto, r), r > R, at initial time uq. 
(We take uq = 0). At the mirror, we set <&(u, R) = A = constant, with the gauge condition 
that oo) = 0. We adopt the coordinate condition (3(u,R) = 0. The condition that the 
metric match continuously to a flat interior for r < R requires V(u, R) = R. (The normal 
derivatives of $ and the metric do not match continuously across r = R in accord with the 
sheet stresses associated with the mirror). 

With these conditions the scalar field and metric components have a unique future evo- 
lution. The resulting metric does not have an asymptotic Minkowski form at X + . This 
is characterized by the quantity H(u) = (3(u, oo) which relates Bondi time t at J + to the 
proper time r = u at the reflecting boundary according to dt/dr = e 2H . Bondi time is 
the physically relevant time for distant observers. The Bondi mass of the system can be 
expressed in either an asymptotic or integral form [6]: 

M(u) = \e-™r\Y-), 

poo 

= 2vr / e 2 ^r 2 (^, r ) 2 dr. (5) 

JR 

(Note that the mass content of a null hypersurface only depends upon field values in the 
hypersurface and is independent of their time derivative). 

In the flat space limit, (2) and (3) imply (5 — and V = r and (4) reduces to the wave 
equation in null-radial coordinates: 

U^g := 2g )Ur - g >rr = 0, (6) 

where g = r$. The general solution satisfying the boundary condition for a kink potential 
A is 

g(u, r) = AR + f(u/2 + r) - f{u/2 + R). (7) 

The first term describes a static monopole; the second, an incoming wave; and the third, 
its reflection off the mirror. At late retarded time, the solution asymptotically approaches 
the static equilibrium state $ — > AR/r. This state is also the solution of the variational 
problem for the energy, SM = subject to the constraint that 5$ vanishes at the mirror 
and at infinity. In this flat space example, 

poo 

M = 2n (§ >r ) 2 r 2 dr (8) 

JR 

and its variation is given by 

poo 

SM = -4tt / [$ r r 2 ] >r 5$dr. (9) 

JR 

In the flat space case, the static solution $ = A/r is the unique extremum satisfying the 
boundary conditions (which rule out the trivial solution $ = constant). Furthermore, the 
second variation of the mass about this state is positive, so that it represents a stable 
equilibrium, which the system approaches after radiating its excess energy to infinity. 
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In Sec. II, we show that for small kink potentials these features extend to the curved space, 
self-gravitating case. But, above the critical kink potential, equilibria do not exist. We show 
that this critical potential marks a turning point instability about which the equilibrium 
configurations bifurcate into stable and unstable form. In Sec. Ill, we present numerical 
simulations of the evolution of stable and unstable kinks. 



II. ENERGETICS OF SELF-GRAVITATING KINKS 

At the initial retarded time, the null data $(0, r) and the constrained kink potential 
A = R) uniquely determine the future evolution of the system. In order to discuss the 
energetics of these states it is convenient to describe the corresponding configuration space 
in terms of $(V) = $(0,r(V)), where V is considered the independent variable. Here V, as 
determined by (2) and (3), is the affine parameter along the outgoing null cone. The affinc 
freedom is fixed by the boundary conditions V{R) = r(R) = R and r'(R) = 1, where we 
denote dyf = f. Then (2) and (3) can be reexpressed as r = /, where / satisfies 

/" = -4tt($') 2 /, (10) 

which determines r(V) uniquely in terms of the initial data. The mass of of the state $(V) 
can be obtained by inserting the solution r(V) of (10) into the first version of (5), which 
gives 



M = -(r-Vr') 
2 V 1 



(11) 

when reexpressed with V as the independent variable. 



A. Static equilibria 

We now show that the asymptotically flat static solution \1> of the Einstein-Klein-Gordon 
system is an extremum of the energy, subject to a fixed kink potential. This solution [7], 
which is the analog of the 1/r solution in a Minkowski background, can be obtained in null 
coordinates by setting $ u = in the wave equation (4). This gives 

rVtf >r = const, (12) 

whose solution, after using (2) and (3) to eliminate the r-dependence, is 

«^ = v^ <+*£Vl> < 13 > 

with r(V) = 7> given by 

r% = e - 4atanha [V + R(e- 2a - l)] 1 - tanha [V + J R(e 2Q - i)]^***". (14) 
Here the integration constant a determines the kink potential, 
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M*) = r a r ■ (15) 
'7r cosh a 



The space-time has a naked singularity when analytically extended to r = [7]. The Bondi 
mass of this solution is 

M*(a) = 2 J Rsinh 2 ae- 2atanha (16) 

In order to investigate the behavior of the Bondi mass with respect to variations about 
\P, we consider the solutions of the second order differential equation (10) for an arbitrary 
configuration $. The solution r is uniquely determined by the boundary conditions r(R) = R 
and r'(R) = 1 and has asymptotic behavior r ~ e~ 2H V + 2M + 0(1/V) as V — > oo. We 
define a second independent solution r by requiring asymptotic behavior r ~ V + 0(1/V). 
For a general configuration this normalizes the Wronskian of these solutions to the Bondi 
mass, 

2M = rr' - rr'. (17) 

Now, for a fixed kink potential, consider the e-dependent family of perturbed configura- 
tions $ = ^ + e(f), subject to the boundary condition (f)(R) = 0(oo) = 0. For e = 0, we have 
r = ?> and r = r^, with given by (14) and 

r 2 = [V + R(e~ 2a - l)] 1 + tanhQ [\/ + J R(e 2Q - i)] 1 -*^. (is) 

To first order in e, 5r = r — satisfies the perturbed version of (10): 

Sr" = -An5r{^') 2 - Mer^'cj)'. (19) 

This linear differential equation for 5r is identical to (10) except for a inhomogeneous term. 
Its solution can therefore be expressed in terms of the two independent homogeneous solu- 
tions and t# as [8] 

47TC f V 

5r = -— / dW[ry(W)Ty(V) - T^(W)r^(V)}r^(W)m'{W)<P'{W). (20) 

Mq, J R 

Note that this satisfies the proper boundary conditions 8r(R) = 8r'(R) = 0. By inserting 
this solution for Sr into (11), we obtain for the first variation of the Bondi mass 



M-Mn 

IR 



poo 

= 4ire / dVT*r^'(j)'. (21) 

JR 

But direct calculation gives 

r*r^' = _ ^ sm ^V 2 Q tan ha = constant (22) 
V 71 " 

so that (21) integrates to yield M - = 0(e 2 ). 

Thus the static solutions are equilibrium configurations. For sufficiently small values of 
a, the second variation of the Bondi mass will be positive, as in the flat space case. In that 
regime, small perturbations can be dissipated by mass loss due to scalar radiation. However, 
if there exists a perturbation which lowers the Bondi mass the static equilibria cannot be 
stable since further radiative mass loss would drive the system away from equilibrium. We 
next see that this indeed occurs in the strongly nonlinear regime. 
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B. The turning point instability 



The one-parameter family of static equilibria has kink potential Ay (a), given by (15), 
which increases monotonically with a from Ay(0) = until it reaches a maximum at the 
turning point a c ~ 1.199 satisfying a c tanho; c = 1. Above a c , Ay (a) monotonically de- 
creases to as a — > oo. Thus, below Ay(a c ), there are two static equilibria for each kink 
amplitude. Similarly, the mass My (a) = increases monotonically from My(0) = to a 
maximum at the same turning point a c and then decreases monotonically to the black hole 
limit, My (a) — > R/2, as a — > oo. 

The quantity k, defined by dMy/da = ndAy/da, plays an important role in the stability 
of this system. Explicitly, 

k = Ay^Rsmh(a)e- 2atanha . (23) 

Denoting / = df /da, it is easy to verify that kAy pa .65R > for a = a c . 

These are precisely the criteria for the application of a theorem [2] regarding the onset 
of instability at a turning point in a one-parameter family of equilibria. We present here 
a finite dimensional version of the argument behind the theorem as applied to the present 
case. Proofs for both the finite dimensional and function space cases are given by Sorkin [2] . 

Consider then the configuration space {$(A X 1 , X n )} whose elements for each choice 
(A, X 1 ) represent a function <$>(V;A, X 1 ) with boundary conditions $(R;A, X 1 ) = A and 
$(oo; A, X 1 ) = 0. Let ty(a) represent the one parameter family of static equilibria, regarding 
(Ay, X 1 ) as functions of a. We now show that for some 5a > all the equilibria in the range 
a c < a < a c + 5a are unstable. 

The proof is based upon the equilibrium conditions, My t A = k and My^ = 0, where 
f,A = 9a f and f\i = dxtf- Differentiating these conditions with respect to a along the 
equilibrium sequence gives 

k = My AA Ay + My Ai X l (24) 

and 

= My :i X l = My M AyX { + A, (25) 

where A = My^X l X^ represents the second variation of the mass with respect to the 
perturbation X 1 . Equation (25) implies that A = at the turning point where Ay — 0, , so 
that this mode has neutral stability. Taking a further derivative of (25) and evaluating at 
the turning point, we obtain 

A = -My M AyX { = -kAy. (26) 

As remarked above, kAy > at the turning point so that (26) implies the instability of this 
mode in a neighborhood a > a c . 
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C. Unstable modes 



Since there are no other turning points, these general considerations suggest that the 
equilibria \& are stable below a c and unstable above a c , although the theorem only strictly 
implies the onset of instability in a neighborhood of a c . In order to provide further insight, 
we now consider specific behavior of the second variation of M^. 

For an arbitrary perturbation, $ = \l/ + e0, we can follow the approach in Sec. II A which 
led to the expression (21) for the first variation of the mass and proceed further to determine 
the second variation. At order e 2 , this leads to 

[•CO 

M-M* = 2ne 2 / dW*(r*0'0' + 25rV'(f)'), (27) 
Jr 

where 5r is given by the integral (20). Because Sr depends nonlocally on the perturbation 
0, (27) does not give straightforward information regarding the sign of the second variation. 
However, the first term in the integrand is positive definite and depends locally on the 
perturbation. Therefore, for perturbations confined to a sufficiently compact region, the 
first term dominates the Sr term and the second variation is positive. This leads us to 
consider large length scale perturbations in the search for unstable modes. 

One such perturbation can be obtained from the static solution by combining a change 
in a with a translation, 

5<f)(V) = V(V + 5R;a + 5a) - V(V;a), (28) 

with 5a and SR adjusted to set S<p(R) = so that the kink potential is held fixed. This 
requires 

p( p — 2(a+Sa) „4acosh(a+<5a)/ cosh(o) p 2(a+5a)\ 

^ = — — A h , ^ w ,< . -■ (29) 

For this perturbation, the change in M can be worked out analytically to second order 
(using Maple) from (17), although the final expression is too lengthy to present here. The 
important point is that the second order variation of M is positive for a < a c , vanishes 
precisely for a = a c and turns negative in a neighborhood a > a c . This corroborates our 
expectation that perturbations are stable below a c and go unstable as we pass through the 
turning point. (It is curious that beyond the turning point M — again becomes positive 
for values of a greater than approximately 2.0). 



III. GRAVITATIONAL COLLAPSE OF UNSTABLE KINKS 

We now examine the behavior of this system numerically, using a null cone evolution 
algorithm for nonlinear scalar waves developed in Ref's [9] and [4]. The algorithm is based 
upon the compactified radial coordinate x = r/(R + r), so that X + is represented by a finite 
grid boundary, with x — 1/2 at the mirror and x — 1 at X + . The code has been tested to 
be globally second order accurate, i.e. the error in global quantities such as the Bondi mass 
is 0(Ax 2 ) in terms of the grid spacing Ax. 
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In the case of a mirror boundary at zero potential (no kink), as studied in Ref. [3], a 
weak scalar field radiates completely to infinity, so that the final mass of the system is zero. 
However, above a critical strength, the scalar field undergoes gravitational collapse to form 
a horizon. In this case, some of the scalar energy is radiated to infinity and the remainder 
crosses the horizon and contributes to the final black hole mass. The mirror itself must fall 
into the horizon for otherwise it would continue to reflect the scalar field until all scalar 
energy were radiated to infinity. Near the critical strength, the sensitivity of the final mass 
is somewhat analogous to the critical behavior studied by Choptuik [10,11], except there is 
now a mass gap because the final black hole must have a mass larger than R/2 in order to 
contain the mirror. In the present case of a mirror with nonvanishing kink potential, we 
would expect this behavior to be modified in several ways. 

First, we consider kink potentials A < A^(a c ). In this case, two static equilibrium states 
ty(a-) and are possible, a_ < a c < a + , with the former expected to be stable and 

the latter unstable. The double valued behavior of M(A) is graphed in Figure 1. It shows 
that Mq,(a + ) > M^(aJ). Consequently, the following scenarios are expected. A state close 
to ^(cl) should return to ^>(a ) but a state close to ^(a+) could either evolve toward 
^(a) or collapse to a black hole. These scenarios have been confirmed by the following 
numerical simulations. (All these simulations are run setting R — 1). 

With kink potential A = A^,(aJ) ~ 0.36563 for «_ = 1, we study the stable equilibrium 
by evolving initial data of the form $ = ^(a_) + A(l — 2x)/2r. For values of A less than a 
critical value A c ~ 0.1929, the system returns to the equilibrium state \1/(q;_) at late times. 
Fig. 2 graphs the corresponding time behavior of the Bondi mass showing its asymptotic 
approach to the equilibrium value M — > M*(a_) .60220. For A > A c , the system collapses 
to form a black hole. 

Next, for the same kink potential A 0.36563, we study the evolution of the unstable 
equilibrium by evolving initial data of the form <£> = ty(a + ) + A(l — 2x)/2r. (Here a + ~ 
1.42429). Fig. 3 graphs the corresponding behavior of the Bondi mass. In this case, because 
of the instability of the static equilibrium ^(a + ), the system undergoes a dynamical change, 
even when A = because of the perturbation introduced by discretization error in the initial 
data. For small A, the state ^(a + ) makes a transition to the stable state \l/(a_). But again 
there is a critical value A c ~ (1.311) 10~ 4 above which the final state is a black hole. 

Now consider a kink potential A > A c , for which no static equilibria exist. With these 
boundary conditions, we would expect any initial state to undergo collapse to a black hole. 
We explore this by considering initial data of the form $ = 2(A C + A)/(r + R). Figs. 4 
and 5 show the resulting behavior of g as a function of x at several times, for A = 0.3 and 
A = 10~ 2 , respectively. At late times, these states indeed form a black hole, as evidenced by 
the flattening of g, outside the cusp being formed at the black hole radius, as the exterior 
field sheds its "hair". 

In summary, the behavior of this system is in complete accord with theoretical expecta- 
tions. 
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FIGURES 
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FIG. 1. Bondi mass as a function of amplitude, for the static equilibria. 
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FIG. 2. Bondi mass as a function of Bondi time for perturbations of a stable equilibrium. 
Different curves represent the different values of the perturbation parameter A: 0.15 (continuous 
line); 0.18 (dashed line); 0.19 (small-dashed line); 0.1925 (dotted line); 0.1928 (dot-dashed line). 
The system collapses to a black hole for A greater than A c srf 0.1929. 
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FIG. 3. Bondi mass as a function of Bondi time for perturbations of an unstable equilibrium. 
The different values of the perturbation parameter A are: (continuous line); 10 -4 (dashed line); 
(1.2) 10~ 4 (small-dashed line); (1.3) 1(T 4 (dotted line); (1.31) 1(T 4 (dot-dashed line). For these 
values the system makes a transition to the stable equilibrium. The system collapses to a black 
hole for A c « (1.311) 10" 4 . 
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FIG. 4. Evolution of g. The initial data is g(0,x) = 2x(A c + A), for A = 0.3. The upper line 
represents u = 0; and the lower one, u = 1.19. The system collapses to a black hole. 
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FIG. 5. Evolution of g. The initial data is g(0,x) = 2x(A c + A), for A = 10 2 . The upper line 
represents u = 0, and the lower one u = 15.33. The system collapses to a black hole. 
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